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ABSTRACT 
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^ , Early results from the Herschel Space Observatory revealed the water cation HtO^ to be an abundant ingredient of the interstellar medium. Here 

\^ we present new observations of the HtO and H2O+ lines at 1113.3 and 1115.2 GHz using the Herschel Space Observatory toward a sample of 

high-mass star-forming regions to observationally study the relation between H2O and H^O^. Nine out of ten sources show absorption from 

H2O* in a range of environments: the molecular clumps surrounding the forming and newly formed massive stars, bright high-velocity outflows 

associated with the massive protostars, and unrelated low-density clouds along the line of sight. Column densities per velocity component of H2O* 

Q are found in the range of 10'^ to a few 10'^ cm"^. The highest N{}\iO*) column densities are found in the outflows of the sources. The ratios of 

5— ( ' YiiOfHiO are determined in a range from 0.01 to a few and are found to differ strongly between the observed environments with much lower 

(Vj ■ ratios in the massive (proto)cluster envelopes (0.01-0.1) than in outflows and diffuse clouds. Remarkably, even for source components detected in 

j^ H2O in emission, HtO^ is still seen in absorption. 

Key words. ISM: molecules - ISM: clouds - Submillimeter: ISM - Stars: formation 

>\ 

^^ ' 1. Introduction a wide range of evolutionary stages is observed in a variety of 

L-^ ' water lines. One of the H20^ ortho ground-state doublet lines 

nn ■ One of the unique ESA Herschel Space Observatory (Pilbratt jj^^ ^j^.g j„ frequency to the H,0 para ground-state line and 

^. et al. 2010) science fields is the observation of thermal j^ observed as well. This allows us to present here a detailed 

C-^ lines of interstellar water and other hydrides. Hydrides have comparison of water and ionized water column densities in a 

O . small reduced masses, so their rotational lines lie at short j^g^^ ^^^^^^ ^f ^o^^ces to study relative abundance variations 

O . submillimeter wavelengths, which are almost unobservable from ^f H2O+ and H2O in different interstellar environments. 
'""' ' the ground (iPhillips & Vastelll2003l) . Early results from the first 

^ \ months of observations show the scientific potential of these 

k> ^^'^'^^^^ (e.g. Ivan der Tak et ai.|[2oI5 for water in a massive star- 2. Observations and data reduction 

/\ forming region). Interestingly, these early results also revealed 

H the water cation HaO^, which was seen by Herschel for the The sources were observed with the Heterodyne Instrument for 

. - - first time, as an abundant ingredient of the interstellar medium the Far-Infrared (HIFI, De Graauw et al. 2010) onboard the 

(Ossenkopfetal. 2010; Gerinetal. 2010). The ortho ground- i/e«c/ie/SpaceObservatory (Pilbratt etal. 2010) on 2010 March 

state line of HaO^ was even detected in external galaxies and 3-5 and April 17. Double-beam switch observations (throw of 

found to be stro i iger than the para grou nd-state water line 2.5 arcminutes) have been performed in the double sideband 

dWeiss et al.l2010tlVan der Werf et al.ll2010h . These earlv results mode using the 4b receiver band. The pointing coordinates 

indicate that H20^ originates mainly from low-density gas of of the observed sample are given in Table [T] Data were 

diffuse interstellar clouds. taken in two polarizations with the acousto-optical wide-band 

Within the "Water in star-form ing regions with Herschel spectrometer (WBS), which covers 4-8 GHz in four sub-bands, 

(WISH)" (IVan Dishoeck et al.l I20T0I) Herschel key program, a each approximately 1.1 GHz wide. Its Nyquist resolution is 

sample of about 20 massive star-forming regions (SFR) covering approximately 1.1 MHz (0.30 km s"). Four species have been 

observed simultaneously with the WBS: P-H2O (1113.3 GHz, 

* Herschel is an ESA space observatory with science instruments USB) and p-Ha O In - Ooo (1101.7 GHz, LSB), CO 10-9 

provided by European-led Principal Investigator consortia and with (1101.3 GHz, LSB), and o-H20^ In - Ooo (1115.2 GHz for the 

important participation from NASA. strongest HE component, USB). 
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Table 1. Sources observed in the 1110 GHz setup. 



G34.26+0.15 
1 . 1 1 45 1 0^ 1.114 10° 1 . 1 1 35 1 0° 



Source Ra (J2000) Dec Vi,,- Lboi 

(hms) (°"') (kms-') (Lq/IO'*) 



AFGL2591 

DR21(0H) 

G29.96-0.02 

G31.41+0.31 

G34.26+0.15 

NGC7538-IRS1 

W3-IRS5 

W33A 

W43-MM1 

W51NE1 



20 29 24.7 
20 39 00.8 
18 46 03.8 
18 47 34.3 
18 53 18.6 
23 13 45.3 
02 25 40.6 
18 14 39.1 

18 47 47.0 

19 23 40.0 



+40 11 19 
+42 22 48 
-02 39 22 
-01 12 46 
+01 14 58 
+61 28 10 
+62 05 51 
-17 52 07 
-01 54 28 
+ 14 3051 



-5.5 
-4.5 
+98.7 
+98.8 
+57.2 
-57.4 
-38.4 
+37.5 
+98.8 
+59.0 



2.0 
1.7 
12 
18 
28 
20 
17 
1.0 
2.3 
10-100 



The system temperatures for our data were around 350 K. 
Integration time was 601s. Calibration of the raw data onto Ta 
scale was performed by the in-orbit system (Roelfsema et al 
2010); the conversion to Tmb was done with a beam efficiency of 
0.7. The Herschel full-beam-at-half-maximum at this frequency 
was assumed to be the theoretical one (20"). Currently, the flux 
scale is accurate to 5%. An rms of 90 mK has been reached. 

Data calibration was performed in the Herschel Interactive 
Processing Environment (HIPE) version 2.8. The velocity 
uncertainty in the current version of the pipeline is up to 2 
km s ', depending on target direction and observation epoch. 
Further analysis was done within the CLASS package. After 
inspection, the data from the two polarizations were averaged 
together. The continuum level in the data was divided by 2, 
because the original calibration was done for the line emission 
originating from only one receiver sideband. 



3. Results 

In Fig. [U the DSB WBS spectrum towards G34.26+0.15 is 
shown as a typical example for the data that are analyzed in this 
study. H2O+ is detected in aU sources except NGC7538 IRSl 
and in all cases seen in absorption, similar to th e previous 
detections dOssenkopf et alJl2010l: lOerin et al.ll201(l ). H2O and 
Hj^O on the other hand show in many cases both absorption 
and emission line components. Many sources show several 
velocity components owing to absorption from diffuse clouds 
from different spiral arms on the Unes of sight (LOS), which 
complicates the interpretation of the H20^ spectra because 
of its complex hyperfin e structure (HFS, see [Strahan et id] 
ri986t lMurtzetalTll998l for details). Several sources show 
saturated H2O absorption down to the K level (e.g. Figs. [1} 
O, demonstrating that the continuum level in the spectra is 
measured reliably and that the sideband gain ratio is 1 . This is 
especially important because most of the analysis is based on 
the detected absorption features for which the line-to-continuum 
ratio is the relevant observing parameter. Most of the sources 
show broad wing emission (AV > 10 km s"') as an indication 
of powerful outflows, in H2O mostly seen in emission, while 
in H20^ outflows are detected as broad (due to the blended 
HFS) blueshifted absorption features in front of the strong dust 
continuum emission. 



4. Analysis 

Because all lines for a given source are observed simultaneously 
in the DSB spectra, they will only have a small relative error 
in their intensities and velocities independent of the calibration. 




1.1015 10 1.102 10 

Rest Frequency (MHz) 



Fig.l. Example DSB spectrum at 1110 GHz towards G34.26 
showing all lines covered in this spectral setup at the Vlsr of the 
source. The redshifted H2O and H20^ features are caused by the 
line-of-sight absorption components. The lower and upper scales 
give the LSB and USB frequency scales, respectively. 



Here, the rest frequencies from iMiirtz et al] (1 19981) are used, 
which have a quoted accuracy of 2 MHz. A comparison of the 
H2O and H20^ velocities from LOS absorption features shows 
a small scatter of +2 km s"', hence the resulting accuracy of the 
measured H20^ velocities is sufficient to associate them with 
known velocity components of the observed sources. 

To separate the various and often blended velocity 
components and the HgO^ HFS, we used the XCLASS 
fitting tool (Comito et al. 2005, and references therein), which 
allows us to obtain multi-component LTE fits of emission 
and absorption component s and which takes the HFS - which 
extends over 40 km s ' (iMiirtz et al.l Il998h - into account. 



■s 2.5 




50 
Velocity (km/s) 

Fig. 2. XCLASS fits to the water (lower panel) and water ion 
line profiles (upper panel) in W43 MMl shown in red solid lines. 
The fit to the Hj^O line (blue spectrum) is shown in green. An 
H20'^ fit without a broad outflow component is shown with red 
dots. The systemic velocity of the source is indicated by a dotted 
line. The strong emission line in the H20^ 
(10-9) from the other sideband. 



spectrum is CO 
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W51N El 




Table 2. H2O, Hj^O and H20^ fit results of velocity components 
associated with the massive star-forming clumps W43 and W5 1 
shown in Figs.|2]and[3j 



50 
Velocity (km/s) 



100 



Fig.3. SameasFig.|2]forW51NEl. 



The input parameters for the fits are the excitation temperature, 
column density, source size, source velocity and velocity 
width. For the absorption components, a background brightness 
temperature of 15-25 K was derived from the measured 
continuum temperatures assuming source sizes as given in 
Table 12 

The H2O absorption components with high velocity offsets, 
which are likely diffiise LOS clouds, were fitted with an 
excitation temperature of 2.7 K. However, the diff'use Galactic 
background radiation might increase the excitation temperatures 
of the submillimeter lines to values of about 5 K. For absorption 
components that are likely associated with the massive SFRs 
we used a fixed value of 5 K, which is clearly below the 
background temperature, to get an absorption in the fit. However, 
the corresponding column densities for temperatures below 10 K 
depend only weakly on the assumed Tex- For the emission 
components we used a fixed value of 50 K. The source size 
is assumed to be much bigger than the beam. Only for the 
57.5 km s"' component toward W51N El a size of 35" is 
assumed. Judging from its strong Hj^O absorption, it has a very 
high H2O optical depth, but still does not absorb the continuum 
down to K and therefore requires a filing factor smaller than 
one. For the column density calculations we used an ortho-to- 
para ratio of 3: 1 for H20^. 

The fit results obtained for H2O were used as the starting 
point to fit the H20^ spectra, which are more complex owing 
to the H20^ HFS. because no emission is seen in H20^, 
we changed any component seen in emission in H2O into 
an absorption component for H20^ by lowering its excitation 
temperature to 5 K. The physical parameters for each component 
were then fine-tuned to fit the observed H20^ spectra. In cases 
where the corresponding component was not detected in H20^, 
the highest column density consistent with a non-detection was 
chosen to derive an upper limit. In a few cases, the observed 
H2O components were not sufficient to account for all the H20^ 
absorption. In Fig. |2] e.g. there is no indication of an H2O 
outflow component and the blue-shifted LOS absorptions are 
quite narrow, while the H20^ absorption is very broad, even 
considering its HFS, so that an additional blueshifted, broad 
"outflow" component (AV > 10 km s"' ) was added to reproduce 
the observed H20^ spectrum. 



Mol. / Source 


Tev 


A'/IO'- 


AV 


V,.,- 




(K) 


(cm-2) 


(Ions-') 


(kms-') 


H2O/W43MMI 


2.7 


6.0 


2.0 


7.0 




2.7 


9.0 


2.0 


10.5 




2.7 


5.0 


2.0 


13.0 




2.7 


3.0 


2.0 


66.5 




2.7 


2.0 


2.0 


70.0 




2.7 


7.5 


2.0 


78.0 




2.7 


3.5 


2.0 


81.0 




5.0 


3.0 


20.0 


87.0 




5.0 


200.0 


4.0 


97.0 


H'^O 


5.0 


5.0 


6.0 


99.0 


H2O+ 


2.7 


3.0 


4.0 


8.0 




2.7 


1.2 


4.0 


11.5 




2.7 


1.5 


4.0 


14.0 




2.7 


3.0 


2.0 


64.5 




2.7 


2.0 


2.0 


68.0 




2.7 


1.0 


2.0 


78.0 




2.7 


3.0 


2.0 


80.5 




5.0 


6.0 


4.0 


96.0 




5.0 


15.0 


20.0 


87.0 


H2O/W5INEI 


2.7 


6.0 


3.0 


6.0 




2.7 


1.0 


5.0 


11.0 




2.7 


0.3 


2.0 


24.5 




2.7 


3.5 


1.0 


45.0 




2.7 


1.5 


2.0 


48.0 




5.0 


400.0 


5.0 


57.5 




50.0 


100.0 


40.0 


59.5 




2.7 


250.0 


7.0 


64.0 


H^**0 


5.0 


3.0 


5.0 


57.5 




2.7 


1.5 


7.0 


64.0 


H2O+ 


2.7 


6.0 


5.0 


6.0 




2.7 


2.0 


5.0 


11.0 




2.7 


0.5 


2.0 


22.5 




2.7 


0.3 


1.0 


45.0 




2.7 


1.0 


2.0 


48.0 




5.0 


17.0 


20.0 


50.0 




5.0 


0.3 


5.0 


57.5 




2.7 


4.0 


4.0 


64.0 



Two examples for the resulting fits are shown in Figs.|2]and 
[3] The corresponding fit parameters for these sources are given 
in Table|2] The range in A^(H20+) of 10'^ to a few 10'^ cm^^ is 
much smaller than the range in A^(H20) of a few 10'^ to several 
10'^ cm"^. Similar results are obtained for other sources in our 
sample and will be discussed elsewhere and in the following 
section. 



5. Discussion and conclusions 

An overview of the fit results is shown in Fig.|4] in which we plot 
for each velocity component the corresponding H2O and H20^ 
column densities. Some of the lower column densities are upper 
limits while some of the high H2O column densities are lower 
limits because of the saturation in the H2O absorption lines. In 
these cases the true column density, estimated from the observed 
Hj^^O lines and assuming a "'0/"^0 ratio in the range from 250- 
560 (Wil son & Rood 1994) wifl be higher by up to a factor 5. 
Figure |4] shows that compared to the properties of d i ffuse clouds 
(see also iGerin et all I2OI0I: lOssenkopf et"an I20IOI: IWeiss et al.l 
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Fig. 4. Comparison of H2O and H20^ column densities of 
different line components in the spectra: the black squares 
give diffuse lines of sight, red triangles the envelopes of the 
massive star-forming clumps, and green filled dots the outflow 
components (AV > 10 km s '). As a reference, the solid 
Une shows N(H20^)=N(H20). Some of the low (high) column 
densities are upper (lower) limits and are indicated by arrows. 



120 lOl) . in which A^(H20+) is closer to A^(H20), significantly 
higher A^(H20)/A^(H20^) is found for velocity components 
likely originating from the envelopes and bright high-velocity 
outflows of the massive star-forming clumps. Even in the 
outflows and envelopes of the massive star- forming regions there 
is still a variation of A^(H20^). The highest A^(H20^) column 
densities are found in the outflows. 

After t he first indication that H20^ is present in the outflow 
of DR21 JOssenkopf et al.l l20 1 Ol) . it is clearly detected here in 
several outflows. Some outflows are even more prominent in 
H20^ than in H2O (e.g. outflow components above the diffuse 
cloud region of Fig. |4|i, although in some cases it is difficult to 
disentangle the blueshifted outflow and spiral arm absorption. 

In Fig. |4]one cloud treated in the litera ture as diffuse LOS 
cloud toward W51 (64 km s ' component. ISolUns et al.l 120041) 
has a high H20/H20^ ratio of 50, which indicates that it might 
rather be associated with the W5 1 clump envelope. Strong ortho- 
wa ter absorption o f this c omponent was already seen with SWAS 
by iNeufeld et all (l2002h . Hence the H2O/H2O+ ratio might be 
useful for the classification of absorption components. 

A key result from our observation of a quite significant 
sample is that H20^ is always seen in absorption, even when 
originating from outflows or envelopes, in which water is 
seen in several cases in emission. This finding might point 
to the origin of H20^ from low-density environments, but 
alternatively, because H20^ is a highly reactive molecular ion 
that reacts rapidly with H2 and electrons, inelastic collisions are 
very ineffective at exciting it into emission, regardless of the 
density. 

The new detections of H20^ toward Galactic star-forming 
regions presented here together with the recent detection of 
OH^ from the ground (Wvrowski et al. 2010 ) and from space 
dGerin et al.1 l2010t iBenz et al., ,2010: ,Bruderer et al.1 12010) are 



an important confirmation of the gas-phase route to water. 
The H20^ lines are stronger than the H30'^ lines in the same 
sources (some of is caused by differences in spectroscopic 
properties), which is surprising, because H20^ is expected to 
react fast with H2 into H^O^, which recombines with electrons 
to produce H2O. This puzzle is even more pro nounced in 
recent Herschel/HIFI observations in diffuse clouds (Geri net aP 
2010) and active galactic nuclei (W eiss et al. 2010), where 
H20^ is even more abundant than H2O itself. While for 
AGN strong X-ray and/or UV radiatio n is likely to dominate 
the chemistry (IVan der Werf et al]|2010h . environments without 
radiation sources require other solutions. One solution might be 
that in the outer envelopes of the massive SFRs H2O is freezing 
out onto grains whereas -in the case of positively charged grains 
- H20^ is much less affected by freeze-out. 

The main destruction routes of H20^ are dissociative 
recombination (into OH) and reaction with H2 (into HsO^ 
and H2O). The high H20^/H20 ratio observed in the diffuse 
components implies that the first channel is faster than the 
second. In gas where all hydrogen is in molecular form, 
the electron fraction is lO"'' at most when all carbon is 
ionized, which is not enough to make recombination faster 
than the reaction with H2. Our observations therefore imply 
that a significant fraction of the hydrogen in the outflows 
is in atomic form. The same conclusion applies to diffuse 
clouds, where UV radiation causes partial dissociation of 
H2 (Gerin et al. 2010), an d also to AGN, where X-rays are 
responsible dVan der Werf et al.ll2010l) . For molecular outflows 
the most likely mechanism to dissociate H2 is by fast (J-type) 
shocks. The required shock velocities of 30^0 km s ' are easily 
reached in the powerful outflows of the sources. 

Models of dense PDRs (e.g. [Sternberg & Dalgarnol 1 19951) 
predict OH^ + H2 — > H20^. In the outflow-walls scenario 
(Bruderer et al. 2009) this leads to a thin PDR layer along 
the outflow wall, where FUV heats and ionizes the gas. New 
modeling of Bruderer et al. (2010, subm.) of hydrides including 
H20^ predicts the abundance of H20^ to be enhanced by four 
orders of magnitude along the outflow compared to the envelope, 
which then could explain the high column densities of H20^ in 
the outflow components. 
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